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In this work, we show that it is possible to study the notion of geodesic deviation equation in 
f(T) gravity, in spite of the fact that in teleparallel gravity there is no notion of geodesics, and 
the torsion is responsible for the appearance of gravitational interaction. In this regard, we obtain 
the GR equivalent equations for f(T) gravity which are in the modified gravity form such as f(R) 
gravity. Then, we obtain the GDE within the context of this modified gravity. In this way, the 
obtained geodesic deviation equation will correspond to the f(T) gravity. Eventually, we extend the 
calculations to obtain the modification of Matting relation. 

PACS numbers: 04.50.Kd, 04.20.Jb, 04.20.0,95. 


I. INTRODUCTION 

The fundamental equation of Einstein geometrodynamics and other metric theories of gravity is the Geodesic 
Deviation Equation(GDE) Q}. It connects the spacetime curvature described by the Riemann tensor with a measurable 
physical quantity, namely the relative acceleration between two nearby test particles. This equation describes the 
tendency of free falling particles to approach or recede from one another while moving under the influence of a spatially 
varying gravitational field. Actually the presence of this kind of tidal force will cause the trajectories to bend towards 
or away from each other which produces relative accelerationQ-Q. Moreover, the important Raychaudhuri equation 
and Mattig relation may be obtained by considering GDE equation for timelike and null congruences. 

One extend gravity theory beyond general relativity is teleparallel gravity (TG). The birth of this gravity theory 
refers back to 1928 |4|. At that time Einstein was trying to redefine the unification of gravity and electromagnetism 
by introducing the notion of tetrad (vierbin) field together with the suggestion of absolute parallelism. In this theory 
the metric gp v is not the dynamical object, instead we have a set of tetrad fields e a (x M ), and instead of the well-known 
torsionless Levi-Civita connection of GR theory, we work with a Weitzenbok connection to introduce the covariant 
derivative [E]. Furthermore, the role of curvature scalar in GR is played by torsion scalar T in the teleparallel gravity. 

Although at the background and perturbation levels TG is completely equivalent to General Relativity, f(T) 
gravity has new structural and phenomenological features. Especially, at the cosmological background it has various 
cosmological solutions which are consistent with the observational data |M3. In addition, by taking spherical 
geometry one can consider the spherical solutions for f(T) gravity EMI- According to these features, f(T) gravity 
is assumed as a viable theory both at cosmological and at astrophysical aspects. Thus, the cosmological and spherical 
solutions in /(T) gravity lead to various viable models that support cosmological observations [Hi, [13 along with Solar 
System tests [l2| in which f(T) must be close to the linear form. In jl8H20l | the authors have studied the Noether 
Symmtery approach in the Friedmann-Lemaitre-Robertson-Walker (FLRW) geometry, to construct some viable f(T) 
functional forms [2ljj . Regarding this line of progress in the context of f(T) gravity, it seems there is still room to 
study some motivating and interesting gravitational and cosmological aspects of /(T) gravity which have not yet been 
studied. Following this idea, in this work we consider GDE in the context of f(T) gravity. 

The GDE has been studied in f(R) gravity theory [23], [23}, so it is appealing to study the GDE in the context 
of f(T) gravity, too. However, there are conceptual differences between GR and TG. In GR which is fundamentally 
based on the weak equivalence principle, curvature is used to geometrize the gravitational interaction and the spinless 
particles follow the curvature of spacetime. In other words, the concept of force is replaced by geometry and the particle 
trajectories are determined by geodesics, rather than the force equation. In TG, on the other hand, the torsion is 
responsible for the appearance of gravitational interaction as a real force. Hence, there is no notion of geodesics in TG. 
In spite of this conceptual difference, one can show that the teleparallel description of the gravitational interaction is 
completely equivalent to that of general relativity [24[. Therefore, it is possible to cast the force equation in TG into 
the form of a geodesic equation in GR and obtain the corresponding GDE in TG. 
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In the present work, we use another approach to obtain the GDE in /(T) gravity. In this regard, first we use the 
method introduced in Refs fihil [28l | to obtain the GR equivalent of f(T) gravity where the field equations are in the 
modified gravity form, such as f(R) gravity. Then, we can benefit of the approach followed in j22| to obtain the GDE 
within the context of f(R) gravity. In this way, we actually obtain the GDE within the context of f(T) gravity. 


II. FIELD EQUATIONS IN f(T) GRAVITY 

Instead of using the torsionless Levi-Civita connection in General Relativity, we use the curvatureless Weitzenbock 
connection in Teleparallelism [E], whose non-null torsion T p ' and contorsion K p are defined respectively by 

- f p „ = e p A {d p e? - d v e$) , (1) 

K p „„ = % - = \{T p p v + T v \ - T p pv \ (2) 

where T p pv is the Levi-Civita connection. Moreover, instead of the Ricci scalar R for the Lagrangian density in 
general relativity, the teleparallcl Lagrangian density is described by the torsion scalar T as follows 

T = S p ^T p pv , (3) 

where 

V" = \ ( K ^ P + 5 P T °"a - d P T “^) • (4) 

The modified teleparallel action for f(T) gravity is given by [29] 

S = 2k j dix \ e \f( T ) + J dAx \ e \£ m, (5) 

where |e| = det (e^) = U—q, 87 tG = k and c = 1. Varying the action (ID with respect to the vierbein vector field e A , 
we obtain the equation [30| 

\dp (eS/n f T (T) - e\T p pX S p f T {T) + S/"d p (T) f TT (T) + , (6) 

where a subscript T denotes differentiation with respect to T and 0 A is the matter energy-momentum tensor, mean¬ 
while all indices on the manifold run over 0,1, 2, 3, and e A form the tangent vector on the tangent space over which 
the metric r] A B is defined. 

On the other hand, from the relation between the Weitzenbock connection and the Levi-Civita connection given by 
the equation ([2]), one can write the Riemann tensor for the Levi-Civita connection in the form 


R p 


p\v 


d x r p pu - d v v p pX + T p aX r^ v - r^r; A 

V^ A - Va K^ v + K p av K\ x - K p aX K^ , 


whose associated Ricci tensor can then be written as 


(7) 


Rpv = V V K% P - V P K p pv + K p av K\ p - K p ap K\ v . (8) 

Now, by using K p pu given by the equation ([2]) along with the relations K^ pv '> cr = = S p ( l/a ' > = 0 and considering 

that S p = 2 K p nu = -2 T p nll one has HHUl 

P — —X7P Q _ n \7PT a — ^P a K 

±lj pLU v ^vpp, y[iv v - 1 - pa u p J MV crpv ? 

R = -T-2S7 p T\„, (9) 

and thus obtain 


G 


g 9 p " T = ~ S° P v, k p™ . 


(10) 
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where G M1/ = R pv — (1/2) g pu R is the Einstein tensor. 

Finally, by using the equation m , the field equations for /(T) gravity in terms of GR quantities, namely the 
equation ©■ can be rewritten in the following form [25l - l28l | 

frG pu + 1 (T/ r - /(T)) g pv + B pv f TT {T) = kQ pv , (11) 

where fr = > /tt{T) = B pv = S up a W a T and Q pv is the matter energy-momentum tensor. Now, this 

equation is in the form of field equation in modified gravity, such as f(R) gravity. 


III. GEODESIC DEVIATION EQUATION IN GR 

Here we start with a little discussion about Geodesic Deviation Equation in general relativity. The geometrical 
meaning of the Riemann tensor is best explained by examining the behavior of neighborhood geodesics. Imagine C\ 
and C 2 are two adjacent geodesics with an affine parameter v on 2-surface S (see Fig.l). The vector field V a = 
is the normalized tangent vector of geodesic C\, and 77 “ = ^2— is the deviation vector of these two adjacent geodesics. 
In total we describe these geodesics with i“(y s). 

Starting with £ v rj a = £ r) V a ([V, rj} a = 0) which leads to V v V v ry“ = V v V v V a and using V x X7 r Z a — V Y VxZ" — 
V [x Y] Z a = R a p 1 sZ^X 1 Y 8 in which Y a = 77 “ and X a = Z a = V a , we can obtain the GDE as follows @ 

^ = ( 12 ) 

As an introduction, here we review briefly the results of finding GDE in GR. We take the energy momentum tensor 
in the form of a perfect fluid 

O pl/ = (P + P)u a up + pg a p, (13) 

where p is the energy density and p is the pressure. The trace of energy momentum tensor is 

0 = 3 p — p. (14) 

Now, we are supposed to calculate R and R pv by looking at the standard form of the Einstein field equations in GR 
(with cosmological constant) 

Rfiu 2 ^9 ^ £ ^-9\xv — /^0/iiy- (16) 

The Ricci scalar and Ricci tensor are obtained as follows 


R = k(p — 3 p) + 4A, 

R^v = k(p + p)u a u f 3 + ^[n(p -p) + 2 A ]g pv . 

Considering these expressions and using the following equation^ 

1 R 

RafijS — ~^{9ct^Rs/3 PotsR'yP 4“ 90sR r yot 9ft'yRsa) ~^r{9ct^9SfH 9ct59 n /p') T G a p^s, 

the right hand side of GDE is 


(16) 

(17) 


(18) 


R a ^sV 0 r,W s 


i(«p + A)e+ ^K(p + p)E 2 


(19) 


where e = V a V a and E = —V a u a . This equation is known as Pirani equation [i|. It is worth mentioning that one 
may derive the Pirani equation for every metric for which the Weyl tensor vanishes. A great deal of important results 
from this equation have been obtained including some solutions for spacelike, timelike and null congruences [3lJ . 
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IV. GEODESICS DEVIATION EQUATION IN f(T) GRAVITY 

In this section, we shall obtain the GDE in GR equivalent theory of /(T) gravity. Before going through the details 
of calculations we extract R by taking the trace of m which results in 


R= — [2 Tf(T) - 2 f(T) A f TT S p ^ P T - *©]. 

JT 


Inserting this Ricci scalar into m we obtain 


1 

,M/ ~ h 


\g^{Tf r - / A JttSp p P V p T - «©) - f TT S,p P V p T A 


( 20 ) 


( 21 ) 


By applying the equation (1181) and considering the zero value of the Weyl tensor C aPl s we will find. 
1 


R, 


oi(3 / y5 


where 


^(^ck7©(5/3 9a5®"f/3 “1“ 9/3S®~/a A ( / Av© A U/ 7 ^ A S P pp'SI P Tf r j'T^){.9a'y95fi 9aS9y P ') (22) 


2/t - 


A {9 ( 17 D 5 p 9 adD^p A 9 P 5 D^ a gp^/Dfia'j fj* 


6/t - 


- 2/ + 2 Tf r A f TT S p y p T - «0 (g ai g S p - g aS g 7 p), 


Dp, = —S,ppS7 p Td T . 


After rising a index in the Riemann tensor and contracting with V^g^V 6 we will have 

R a p lS V 0 g^V s = 


(23) 


(24) 


2/t 


K(S"0 S p - SfQtf A gp S &" - ff^0?) + (-/ - k© + Tf T A S p p P S7 p T f TT ){5«g S p - Sfg^p) 


+(6«D sp - SfD^p A gp,D a - g Pl D%)f T V 0 gW s - — - 2/ + 2 Tf T + f TT S p p P S7 p T - kQ 

J o JT L 

{ga^gsp gasg^p^jV 0 g~* U ■ 

Note that Eas. (l22l) and (l24l) hold only if the Weyl tensor vanishes. By considering the equation (fl3l) in the equation 
we find 


R° 


Pit 


1 

2/t 


^{p ~b P) (9a~/'U5'U j P gaS^ 7^/3 A 9p5V"yUa gp^^S^oi) 


(25) 


'2np 2 / 

{g al gs p gaSg^p) ( ^ A Tf^p/S A “>5^ /tT ) “b {ga'yDgp 9adD*y P A 9p5D^ a 9{3jDs a ^/t 


Under the condition of vector field normalization, we have V a V a = e and 

RaprysV^V 6 = [k(p + p)( gait'llpV 13 ) 2 - 2(upV l3 )V( oc u 7 \ + eu a u^) + 

^JT 

i~3~ ~1T~ A ^ MpV p T/ TT - 0 (e ffaT - U a U 7 ) 

H - [(*7a:7-^<5/3 QoiSD'yfi “t~ 9(35^ 7 o; 5^7 ^)/t]^^]- 
Again we rise the first index then contract with g 1 


(26) 
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R a ^sV^r,W 5 = 


By considering E = — V a u a 


—— n(p + p){{upV^) 2 rf‘ - (upV 0 )V a (u 1 ri' 1 ) - {upV^vP^rp) + eu a u 7 g 1 ) + 
2 JT l 

^ ppV^T/rr - 0 (erf - V a (V^)) 

+ [(S°D Sfi - SfD^p + g ps D« - g M D<i)f T }V^V 5 ^ . 
and gaU a = g a V a = 0, the equation (1271) converts into 


(27) 


iT^VVV 5 = 77 ^ K(p + p)E 2 +e 
^}T L 


/ 2np 

{ — 


+ ^% p V^T/ tt - 


L ) 

3 ) 


+ 


2/t - 


\S°D SP - 5<ZD lP + gp S D “ - gp^)f T ]V^V 5 


(28) 


It should be mentioned that all the results obtained in this section do not require the FRW ansatz, rather are valid 
as long as the Weyl tensor vanishes. In the next section, we will use these results to obtain the GDE in GR equivalent 
of f(T) gravity model for FLRW metric whose Weyl tensor is zero. Obviously, our final result will be acceptable 
provided that the GR limit of this model be checked. 


A. GR equivalent method with FLRW background 


We take the standard model line element (FLRW universe) including a(t) and k. respectively as scale factor and 
spatial curvature of the universe, as 


ds 2 = —dt 2 + a 2 {t) 


dr 2 


1 — kr 2 


+ r 2 dd 2 + r 2 sin 2 0d<jr 


(29) 


whose Weyl tensor is zero because of conformal flatness. According to the second section, the expression for the 
torsion scalar for the flat FLRW metric reduces to 


T = -6 ff 2 , (30) 

where H = ^ is the Hubble parameter. Note that the equation (l30l) holds only if a diagonal FRW tetrad is chosen, 
since different tetrads giving the same metric lead to different results in f(T) gravity |32| . Apparently, the torsion 
scalar is thoroughly time-dependant, thus we are supposed to be concerned just about time derivatives of T in D M1/ . 

The vector field normalization implies that V a V a = e and also we have E = —V a u a , g a u a = g a V a = 0, rjou° = 0. 
Not only these mentioned conditions but also the non-vanishing components of S tensor will be used to extract the 
final result for the action of operator on fo. Therefore, after cumbersome calculations we obtain 


0- {[&*D SP - SfD^p + gp S D« - gp 1 Dg)f T ) V^V 5 rf = -24 H 2 Hf TT (E 2 + 2e) V ° 


(31) 


where have used S} 0 = S 2 0 = Sf 0 = —2 H{t). Thus, we can write the reduced expression for R a p 1 sV l3 g 1 V s as follows 


R a p-ysVPgW 6 = -i- ( K (p+p)-24H 2 Hf TT )E 2 + 
2-JT L 


,2 , (^P _ / Tf T \ 

V 3 3 3 /J 


1 , 


which is the generalized Pirani equation. Now, we can write the GDE in f(T) gravity model 


2 „a 


D 2 r( 

Du 2 


1 

2/t - 


(n(p + p)-24H 2 Hf TT )E 2 + e 




V 3 


3 J 


V 


(32) 


(33) 


Apparently, as a result of homogeneity and isotropy of FLRW metric, in this equation only the magnitude of the 
deviation vector 77 “ is changed along the geodesics. Whereas in anisotropic universe, like Bianchi I, we can also infer 
a change in the direction of the deviation vector, as described in |33j. 




















B. Direct method with FLRW background 


Once we assume an FLRW metric (tetrad), then the Riemann tensor can be calculated. This means that the 
geodesic deviation equation m can be straightforwardly written in terms of H and its derivative. We can then use 
the f(T) gravity cosmological equations (PHI) and (PHI) (see below) to connect H and H to the cosmological sources. 
Following this way, we can recover equation (PHI) as follows. By taking radial spatial part, for suitable choice of local 
coordinates, we can write LHS of the equation (PHI) for a = r as follows 

R r p^V p rpV s = R r t-ytV^V 1 + R r ri rV r ^V r + R r e -,eV e rfV e + R r wV+rj W+. (34) 

By considering non-vanishing components of the Riemann tensor for the FLRW metric, we must take 7 = r. Thus, 
the equation (l34l) can be written as 

R r MsV p rfV s = R r + R r rrr + R r 6r e g ee V e V e rf+R r ^ gt+ViV+rf = {-HE 2 +eH 2 ) V r , (35) 

where we have used R 4 R 4 = E 2 , Virf = 0, R r rrr = 0, R r e r e = r 2 a 2 , R r 4 > r <t> = a 2 r 2 sin 2 6 and R r ttr = f- Note that 
similar equations are also obtained for a = 9 and a = <j>. 

According to the equations (PHI) and (PHI) (see below) we have 


H 2 


1 (W_f\ 

2fr\3 6 / ’ 


(36) 


and 


H = -^-(« ( P+p)+ 4 ^T/ tt ) . (37) 

By substituting the equation s (1361) and (1371) into the equation (1351) we can obtain the generalized Pirani equation 
(El by means of the direct approach as 


R a wVtrPV* = 

4/t 


C K{p + p)-2AH 2 Hf TT )E 2 


(2 k P f T/ t n 
+ £ (— ~3 + — ) 


which results in the same geodesic deviation equation (1331) . This means that we have found the same results by two 
different approaches which confirms the validity of the obtained geodesic deviation equation. 


C. Fundamental observers with FLRW background 


Here, we are going to limit ourselves to the fundamental observers. In this particular case, we interpret V a and v 
(the affine parameter) as the four-velocity of the fluid u a and t (the proper time), respectively. Since we are treating 
with temporal geodesics we have e = — 1 and also we fix the vector field normalization by E = 1 which leads to 


R a ^ A v = ± 
2jt 


2 Kp 


f 


—+ Kp — 24 H 2 Hf TT + — 


L 3 


6 J 


(38) 


We know that if rj a = £e a , where e a is parallel propagated along t, then the isotropy results in 


from which we have 


De a 

~Dt 


= 0, 


D 2 ri a d 2 i „ 
Dt 2 ~ dt 2 6 


(39) 


(40) 


By using (fl2l) and (l38l) we can write 











7 


which for the particular case i = a(t) leads to 

a 1 
a fr 


Kp Kp + f 
3 2~ + UH HfTT ~ 12 


(42) 


This equation is nothing but a special case of the generalized Raychaudhuri equation. It is worth to mention here 
that the above generalized Raychaudhuri equation can be obtained by the standard forms of the modified Friedmann 
equations in f(T) gravity model for flat universe [30} 


^ 2 = i ( p+ ^ (2T/t ~ /_t) )’ (43) 

2 H + 3 H 2 =-up- AHTfrr - 2.Hf T + 2 H + Tf T - ^ ^. (44) 

The consistency between the modified Friedmann equations in f(T) gravity for flat universe [30| and the generalized 
Raychaudhuri equation for flat universe (l42l) confirms the fact that the approach followed here is a correct one. 


D. Null vector fields with FLRW background 

In this section, we consider the null past directed vector fields, namely V a = k a ,k a k a = 0, for which the equation 
(El reduces to 


R a wkPrptf = —(*(p + p) - 24 H 2 HfTT)E 2 rt a . 
6JT 


(45) 


Actually, this is Ricci focusing in f(T) gravity as is explained in the following. By considering rj a = rje a , e a e a = 1, 
t a u a = e a k a = 0 and also writing an aligned base parallel propagated = k p X7pe a = 0, we obtain the null GDE 
El new form as follows 


|| = -^(K(p + p) - 2AH 2 Hf TT )E\ (46) 

According to the GR studied in @] all classes of past-directed null geodesics experience focusing if we have n(p+p) > 0. 
Therefore, in a particular case with the equation of state p = —p (cosmological constant) we can not recognize any 
focusing effect. Obviously (l46l) shows the focusing condition for f(T) gravity model provided that 


K{p + p) ^ 24 H 2 Hf TT 

h > Tt 


(47) 


Now, we have an expression (14611 which can be written in terms of the redshift parameter z. To do this, we may 
write 


d dz d 
dv dv dz 


which results in 


(48) 


dv 2 


(dv\~ 2 


1 d 2 v d 

d 2 ' 

\dz) 

\dz) 

dz 2 dz 

+ dz 2 


(49) 


Let us consider the null geodesics for which we have 








.. .ao E dz da 

(1 + z) = — = —- > --=-. 

a E 0 1 + z a 


(50) 


Choosing oq = 1 (the current day value of the scale factor), leads to the following result for the past-directed case 


Thus we obtain 


and so 


dz = (1 + z) — ^dv = (1 + z)—Edv = E 0 H( 1 + z) 2 dv. 
a dv a 


dv 


dz E 0 H(l + z) 2 ’ 


(51) 


(52) 


where 


drv 1 I" 1 dH o 

T 2 = ~ E 0 H(1 + z) 3 [H^ 1 + Z ^~dz + 


dH dv dt dH 1 dH 

dz dz dv dt if (1 + z) dt 


(53) 


(54) 


where use has been made of ^ = E = Eq( 1 + z). From the definition of Hubble parameter we can write 


Using (TT2]) . H becomes 


H = --H z . 
a 


(55) 


H h 


Zl!L-^ + nHH‘f TT -L 

3 2 J 12 


-if 2 , 


(56) 


thus 


d 2 v 

T 2 


E 0 H( 1 + z) 3 
Putting this result in (l49l) . we obtain 


1 + 


-J— (™ + ™- 12 HH>f TT + L 

3 H 2 f T \ 3 2 J 12 


(57) 


d 2 r] 
fa 2 


= E 0 H(l + z) 2 


d 2 ri 
fa 2 ' (1 +z) 


+ 




Finally, using (l46l) the null GDE takes the following form 


dr] 

dz 


d 2 r] 3 
fa 2 + (1 + 2 ) 


1 + 


1 


3 H 2 f T V 3 


up up 


- 12 HH 2 f TT 


12 


drj n(p + p) — 24 H 2 Hf T 

Tz + 


2ii 2 (l + z) 2 / T 


-t] = 0. 


Matter and radiation contributions to p and p can be written respectively as 


K P = 3H 2 n m0 (l + zf + 3H 2 n r0 (l + z) 4 , 


up = if 0 2 U r0 (l + z) 4 , 


(58) 


(59) 


(60) 


where we have used p m = 0 and p r = | p r . By using the equation (l60l) the null GDE equation (l59l) can be written as 


0 + P{H, H,z)^ + Q{H , H, z)r, = 0, 


(61) 


where 
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P(H, H, z) = 


H m o(l + z) 3 + |H r o(l + z) 4 + 12IT 1 (3/t + T g T )[^mo(l + z) 3 + H r o(l + z) 4 + £Ide\ 


Q(H,H,z) = 


/t(i + z) [n m o(i + z ) 3 + ri r o(i + z ) 4 + r^Ds] 

3H m o(l + z ) 3 + 4H r0 (l + z ) 4 + 2 ^ T (H m o(l + z ) 3 + H r g(l + z ) 4 + Qde) 


2/t(1 + z) 2 [Q m o(l + z) 3 + H r o(l + z ) 4 + £Ide] 
in which we have applied the following new form of (1431) 

H 2 = i?Q[f2 m o(l + z) 3 + fij-o(l + z) 4 + fi,Di;], 

where fins has been defined as 


(62) 

(63) 

(64) 


n 


DE = 


H 0 




(65) 


Note that in solving the equation (l6lT) . we must use (13U1) . In order to check for the agreement of the above results 
with those of GR, we take the particular case f(T) = T — 2A. As a result of this choice we have fr= 1 and /tt = 0. 
Furthermore f Ide reduces to 


Mde = 772 
n o 


1 (T- 2A 


T - 2A + T 

6 


A 

3 m 


= Ha- 


which can be written as the Friedmann equation in GR 

H 2 = H 2 [H m ()(l + ~) 3 + H r o(l + z) 4 + Ha]- 
Hence we hnd the reduced expressions P and Q as follows 


|H m o(l + z) 3 + 4H r o(l + z) 4 + 2Ha 
( 1 + z)[S2 m o(l + z) 3 + H r o(l + z) 4 + Ha] 

. . _ 3H m o(l + z) + 4H,.o(l + z) 2 

2[H m o(l + z) 3 + H r0 (l + z) 4 + Ha] 

Eventually, the GDE for null vector fields becomes 

cPrj 7fi m o(l + z ) 3 + 4H r o(l + z) 4 + 2Ha dr/ 3H m o(l + z) + 4H r Q(l + z) 2 
dz 2 (1 + 2)[fi m o(l + z) 3 + H r o(l + z) 4 + Ha] dz 2(fi m o(l + z ) 3 + fi r o(l + z) 4 + Ha) 

In order to obtain the Mattig relation in GR [34), we have to fix Ha = 0, H r o + H m o — 1 which lead to 


( 66 ) 

(67) 

( 68 ) 

(69) 

(70) 


d 2 r] ffi m o(l + z ) 3 + 4H r0 (l + z) 4 dr/ 3fi m o(l + ~) + 4fi r o(l + z) 2 _ . . 

dz 2 (1 + 2)[fi m o(l + z) 3 + H r0 (l + z) 4 ) dz 2(H m o(l + z) 3 + H r o(l + z) 4 )^ 

So, the equation (1BT1) gives us an opportunity to generalize the Mattig relation in /(T) gravity. By considering the 
last result, we can infer the following expression for the observer area distance rg(z) [341 ]: 


r 0 (z) 


/ dA 0 (z) 


V(z')\ z 

dfi 


dr)(z')/d£ | z / =0 


(72) 


where A 0 is the area of the object and also H is the solid angle. Equipped with the d/d£ = i? 0 1 (1 + z) 4 d/dv = 
H( 1 + z)d/dz and setting the deviation at z = 0 to zero, clearly we have 


M*) 


_(Kf)_ 

H(0)dr](z , )/dz'\ z ' =o 


where H( 0) is the evaluated modified Friedmann equation (IffTH) at z = 0. 


(73) 
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V. NUMERICAL SOLUTIONS OF GDE FOR f(T) GRAVITY 

To solve null GDE (lUTl) in f{T ) gravity we should consider /(T) functional forms. One simple assumption in 
simplifying the equation (RTTll is T = 0. However, this reduces our model to Lambda-CDM since Qde = Const , as can 
be seen from the equation (1651) . To find new interesting features of /(T) modified gravity we should consider models 
beyond the T = 0 assumption. Actually, there are variety of f(T) functional forms in the literature, each of which 
having some interesting specific features. Even, some studies have been done to obtain viable /(T) functional forms, 
using the Noether symmetry approach 00. 

Following the above discussion, and for simplicity, first we assume T = 0 which is equivalent to a constant T = 
—6 H 2 = To. Thus, the equation (l62l) can be rewritten as follows 

P + P(H, T 0 , z )¥-+Q(H, To, z)t) = 0, (74) 


where 


P{H, Tq, z) 


+ z ) 3 + |f2 r o(l + z ) 4 + i^hI + 3/:r(Tb)[f2 m o(l + z ) 3 + D r o(l + z) 4 + Qde(To)} 
/t( 2 o )(1 + z)[D m o(l + z ) 3 + Hro(l + z ) 4 + Qde(Tq)] 


(75) 


OtHT 1— 3D m0 (l + z) + 4fi r0 (l + z) 2 

^ ’ °’ z >-2f T (T 0 )[n m0 (i + z) 3 + n r0 (i + z) 4 + n DE (T 0 )Y 

and 


(76) 


O 1 fToMTo) (f(T 0 )+T 0 \\ 

QdM= hI{— --6- ))■ 


To solve the equation (I74f) . one can choose the functional form of the /(Tq) as follows 


/(T 0 ) = aT 0 + /3Tq , 


(77) 


(78) 


where a and f3 are constants. We can solve the equation (1741) numerically which results in rj(z) and r 0 (z) plotted in 
Fig.2, as a function of z. 


I?(z) r i; (z) 




FIG. 2: Plot of the deviation vector magnitude r/(z) (left) and plot of the observer area distance ro(z) (right). The parameter 
values chosen are Ho = 80 km/s/Mpc, fl m o = 0.3, fl r o = 0, Qa = 0.7, and dr/(z)/dz\ z= o = 0.1. 

Now, we assume T ^ 0. Since solving the GDE (1611) is not an easy task for any given f(T) functional form, and 
in order to much benefit of numerical calculations in solving this complicate equation, we just consider the following 
power law form as 
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FIG. 3: Plot of the deviation vector magnitude r/(z) (left) and plot of the observer area distance ro(z) (right). The parameter 
values chosen are Ho = 80 km/s/Mpc, fi m o = 0.3, a = 0.9, b = 0.2, c = 0.7, n = —1 (solid), n = —2 (dotted), and n = 2 
(dashed), fi r o = 0 and dr/(z)/dz\ z =o = 0.1. 






FIG. 4: Plot of the deviation vector magnitude r/(z) (left) and plot of the observer area distance ro(z) (right). The parameter 
values chosen are Ho = 80 km/s/Mpc, fi m o = 0.3, a = 0.9, b = 0.2, c = 0.7, n = 1 (solid), n = 1.2 (dotted), and n = 1.8 
(dashed), fi r o = 0 and dr/(z)/dz\ z= o = 0.1. 


f(T) =aT + bT n + c, 


(79) 


which is well suited in 


most of the cosmological implications of f(T) gravity 111,®, 


or in a slightly modified form as 


f(T) = aT + 


bT n 


1 _|_ 7dn+l) ’ 

where a, 6, c and n are constants. One can rewrite the equation (1441) as 

T k P + 3 H 2 - T f T 


L 

2 


12 H 


4T frr + 2/t 


(80) 


(81) 


Thus, by plugging the equation (IBTTi in the equations (RJ21) and (lf>3l) . we can solve the equation (l6lT) numerically for 
suggested f(T) forms (17^1) and (1501) which results in r](z) and r 0 (z) plotted, respectively in Fig.3 for f(T ) = aT+bT n +c 
and in Fig.4 for /(T) = aT + 1+ jffi„ + i) + c, as a function of 2 . 

For the model T = —6H 2 = To, the behaviour of null geodesic deviation and observer area distance are the same as 
those of ACDM model. This is expected because for T = 0 our f(T) model is reduced to ACDM model. For the other 
suggested f(T) models, namely f(T) = aT + bT n + c, and f(T) = aT + 1+ ^„ +1) + c, the general behaviour of null 
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geodesic deviation and observer area distance are similar to those of ACDM. At small values of z < 0.2, the relative 
deviations and relative observer area distance with respect to the ACDM model are almost ignorable. However, for 
redshifts z > 0.2, the geodesic deviation r/(z) and the area distance r 0 (z) corresponding to the suggested f(T) models 
are rather larger than those of ACDM model. This indicates that the suggested f(T) models predict more acceleration 
than ACDM model, for large values of redshift. 


VI. CONCLUSIONS 

In this paper, we have considered Geodesic Deviation Equation(GDE) in the GR equivalent of f(T) gravity model. 
First, we have calculated the Ricci tensor and the Ricci scalar with the modified field equations in f(T) gravity theory. 
Then, in FLRW universe, the Geodesic Deviation Equation corresponding to these GR equivalent quantities of f(T) 
gravity is obtained. To show the consistency of our approach in constructing the GR equivalent of f(T) gravity, 
the generalized GDE and Pirani equations are recovered for f(T) = T — 2A. We restricted our attention to extract 
the GDE for two special cases, namely the fundamental observers and past directed null vector fields. In these two 
cases we have found the Raychaudhuri equation, the generalized Mattig relation and the diametral angular distance 
differential for f(T) gravity theory. We have also obtained the geodesic deviation rj(z) and the area distance ro(z) 
corresponding to two suggested f(T) models and compared them with those of ACDM model. 
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